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Abstract 

A one-dimensional Hubbard model with nearest and (negative) 
next-nearest neighbour hopping is studied variationally. This allows 
to exclude saturated ferromagnetism for U < U c . The variational 
boundary U c (n) has a minimum at a "critical density" n c and diverges 
for n — > 1. 
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The old question of the existence of ferromagnetism in the single-band 
Hubbard model has found a renewed interest, when rigorous proofs for a 
large-spin ground state were discovered for certain special cases [1, 2, 3, 4]. 
While it is in general very hard to show that the ground state is ferromag- 
netic, it is often rather easy to exclude saturated ferromagnetism, because the 
lowest-energy fully polarized TV-electron state \F) is trivially constructed by 
occupying each of the N lowest single-particle levels by a spin-up electron. To 
exclude saturated ferromagnetism one has then simply to find another state 
of lower energy. This approach has been applied to the Nagaoka problem [5] 
in terms of variational iV-electron wave functions with a single overturned 
spin [6, 7, 8]. We have recently used such an approach to show that for 
non-pathological band structures and dimensions d > 2 the ground state is 
neither fully nor partially polarized in the low-density limit [9]. On the other 
hand a low-density route has been proposed for d = 1 by Muller-Hartmann 
[10]. It is this case that we want to discuss in the present contribution. We 
consider the ti — t 2 Hubbard chain described by the Hamiltonian 

L L 

H = - (*i c !<7 C i+i<r + hc^ci+ia + h.c.) + n n n H (!) 

i = l,CT=t,4 1=1 

For t 2 = we recover the simple Hubbard model with only nearest-neighbour 
hopping, for wich ferromagnetism has been discarded a long time ago by Lieb 
and Mattis [11]. For t 2 ^ particles can pass each other and do not preserve 
a given order. The Lieb-Mattis theorem is then no longer applicable and 
indeed, for U = oo and t 2 < ferromagnetism appears both for one hole in 
a half-filled band [12] and in the limit t 2 — >■ for arbitrary densities (n < 1) 
[13]. These results raise the question about the range of parameters where 
the ground state is ferromagnetic. Before addressing this problem we turn 
to the easier question for which parameter values the fully polarized state is 
unstable. 

The fully polarized state |F) is certainly unstable if the energy is reduced 
by flipping a single spin. We use two different trial states. The first ansatz 
[6], defined by 

\^) = f[(l- V n n n li )cl oi c kFt \F) (2) 
i=i 

removes an up spin from the highest occupied level and puts a down spin to 
the minimum of the single-particle spectrum. It is important to note that this 
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minimum moves away from k = for large negative values of h < — 4 , 

namely to k = ± arccos(-^). It is precisely in this region where according to 
Muller-Hartmann ferromagnetism should appear at low density and infinite 
U [10]. The variational parameter r\ in Eq. (2) reduces the double occupancy 
and steadily increases from to 1 as a function of U. The energy change due 
to the spin flip consists essentially of a gain in band energy and a loss due 
to the Coulomb repulsion [9]. Thus one expects the fully polarized state to 
be unstable for small U. This is confirmed by explicit calculations shown in 
Fig. 1. 

As a second trial state we use the wave function introduced by Edwards 

[14] 

1 L N-1 

ix> = -7=E^n4(oio> ( 3 ) 

V L l = l Q=l 

where 

L 

4t(0 = Y, Va{m-l)c ] m ^ (4) 

m—l 

creates an up-spin electron in an orbital that is determined variationally. The 
variational parameters are the wave vector q and the (N — 1)L amplitudes 
(p a {l)- The ansatz (3) is exact for t 2 = [14]. For t 2 7^ we add the Gutzwiller 
projection operator, which amounts to replace the one-particle orbitals (p a (n) 
by (1 — r)8 n fi)(p a (n). Then the wave function (3) includes the ansatz (2) as 
a special case. The variational energy for orthonormal one-particle orbitals 
reads 

E(q, Wa(l)}) = 2*i cos(g) det + 2t 2 cos(2^) det 5 (2) + 2*i trS (1) 

+ 2t 2 trS( 2 ) + UJ2 ^(0)^(0) (5) 

a/3 

with overlap matrices S^p defined as 

s$ = Eva(0^(i + (6) 

1=1 

We have used the conjugate gradient method [15] to minimize the energy. 
The derivatives have been calculated analytically, and after each iteration 
the orbitals have been orthonormalized using the modified Gram-Schmidt 
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method. As initial orbitals, which have to be close enough to the final wave 
functions ip a (l), we use the N — 1 lowest eigenfunctions for a single-site 
impurity. With this second ansatz, which involves much more variational 
freedom than the first, the instability line is moved to higher values of U, as 
shown in Fig. 1. It is worthwhile to add that a still more refined variational 
wave function [8] yields a similar instability line [16]. 

Our variational calculations delimitate a low U region where the ground 
state is certainly not fully polarized. For the present model and large enough 
negative values of t 2 there is still a large region left where saturated ferromag- 
netism cannot be excluded, as illustrated in Fig. 1 for the particular value 
t 2 = — 0.8 |^i | . In order to find out whether these variational lower bounds for 
the critical value U c are likely to be representative for the exact boundary, 
we have carried out exact diagonalizations for L < 20. The numerical results 
show that for densities around 0.4 and t 2 ~ — 0.5|ti| the exact boundary lies 
typically about 20-50% above the variational bound of the second ansatz, Eq. 
(3). An interesting question is what happens above the "critical density" n c 
where the Fermi surface is reduced from four to two points. Both instability 
lines shown in Fig. 1 have a minimum at this density. In addition, for the 
Edwards ansatz the optimal wave vector q is found to jump from q = ko 
(the location of the minimum in the band structure) for n < n c to q = tt for 
n > n c . In our previous exact diagonalization studies (for an even number 
of particles and periodic boundary conditions) we tentatively associated n c 
with the phase boundary for ferromagnetism. This would indicate that the 
low-density route to ferromagnetism is not directly connected to the Nagaoka 
problem. In the mean-time we have extended our numerical studies both by 
imposing other boundary conditions and by increasing the chain length, using 
the density-matrix renormalization group method. The results obtained so 
far do not reproduce the phase boundary at n c . Instead, at least for U — > oo 
and not too large \t 2 \, ferromagnetism seems to exist for all densities n < 1. 
Further calculations will be needed to establish the full phase diagram. 
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Figure Captions 

Fig. 1. Variational phase diagram of the ti-t 2 Hubbard chain in the n-U 
plane, for t 2 = — 0.8|ti|. Full line : ansatz (2) for L = 1000, dots : ansatz (3) 
for L = 100. 
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